Statistical Ensembles with Fluctuating Extensive Quantities 
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We suggest an extension of the standard concept of statistical ensembles. Namely, we introduce a 
class of ensembles with extensive quantities fluctuating according to an externally given distribution. 
As an example the influence of energy fluctuations on multiplicity fluctuations in limited segments 
of momentum space for a classical ultra-relativistic gas is considered. 
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Successful application of the statistical model to 
hadron production in high energy collisions (see, e.g., 
recent papers [l[ and references therein) has stimu- 
lated investigations of properties of statistical ensem- 
bles of relativistic hadronic gases. Whenever possi- 
ble, one prefers to use the grand canonical ensem- 
ble (GCE) due to its mathematical convenience. The 
canonical ensemble (CE) 0] should be applied when 
the number of carriers of a conserved charges is small 
(of the order of 1), such as strange hadrons [^1, 
antibaryons i, or charmed hadrons i, in other- 
wise large systems. The micro-canonical ensemble 
(MCE) has been used to describe small systems 
(with total number of produced particles less than 
or equal to 10) with additionally fixed energy, and 
momentum, e.g. elementary particle collisions or an- 
nihilation. In these cases, calculations performed in 
different statistical ensembles yield different results. 
Hence, ensembles are not equivalent, and systems are 
'far away' from the thermodynamic limit (TL). 

Measurement of hadron multiplicity distributions 
P{N) in relativistic nucleus-nucleus collisions opens 
another interesting field of investigations. The aver- 
age number of produced hadrons ranges from 10^ to 
10'', and mean multiplicities (of light hadrons) ob- 
tained within GCE, CE, and MCE approach each 
other. One refers here to the thermodynamical equiv- 
alence of statistical ensembles and uses the GCE for 
fitting experimentally measured mean multiplicities. 
However, the number of particles fluctuates event-by- 
event. These fluctuations are usually quantified by the 
ratio of variance to mean value of a multiplicity dis- 
tribution P{N), the scaled variance, and are a subject 
of current experimental activities. In statistical mod- 
els there is a qualitative difference in the properties of 
mean multiplicity and scaled variance of multiplicity 
distributions. It was recently found 0,11, [§| that even 



in the TL corresponding results for the scaled vari- 
ance are different in different ensembles. Hence the 
equivalence of ensembles holds for mean values in the 
TL, but does not extend to fluctuations. 

Statistical mechanics is usually formulated through 
the following steps. Firstly, one fixes the system's ex- 
tensive quantities: volume V, energy E, momentum 
P, and conserved charges {Qi}. Secondly, one pos- 
tulates that all microstates have equal probability of 
being reahzed. This defines the MCE. The CE intro- 
duces temperature T. Each set of microstates with 
fixed energy E (a macrostate) is weighted with the 
Boltzmann factor e~^/^. The probability to find a 
macrostate with energy E in the CE is then propor- 
tional to the number of all microstates with energy 
E times the Boltzmann factor e~^^'^ . To define the 
GCE, one makes a similar construction for conserved 
charges {Qi\ and introduces chemical potentials {/ii}. 
Canonical or grand canonical observables can then be 
obtained by averaging over the CE energy distribution 
or the GCE (joint-) distribution of both energy and 
conserved charges. 

Fluctuations in statistical systems, e.g., multiplic- 
ity distributions P{N) in relativistic gases [zl, H 01, 
are sensitive to conservation laws obeyed by the sys- 
tem, and therefore to fluctuations of extensive quanti- 
ties. For calculation of multiplicity distributions, the 
choice of statistical ensemble is then not a matter of 
convenience, but a physical question. Fluctuations of 
extensive quantities A = (V, E, P, {Qi}) around their 
average values depend not on the system's physical 
properties, but rather on external conditions. One 
can imagine a huge variety of these conditions, thus, 
MCE, CE, GCE, or pressure ensembles [l^ are only 
some special examples. 

A more general statistical ensemble can be defined 
by an externally given distribution of extensive quan- 
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titles, Pa{A). All microstates with a fixed set A 
are taken to be equlprobable. Thus, the probability 
Pmce{N\ A) of finding the system in a macrostate with 
fixed A and additionally fixed multiplicity TV is given 
by the ratio of the number of microstates with fixed 
N and A to the number of microstates with fixed A. 
The construction of multiplicity distributions in such 
an ensemble proceeds in two steps. Firstly, the MCE 
multiplicity distribution, Pmce{N', A), at fixed values 
of the extensive quantities A is calculated. Secondly, 
this result is averaged over the external distribution 
PciA), 



P^{N) = J dAP^iA) P„ 



^iN;A) 



(1) 



The ensemble defined by Eq. ([T]), the a-ensemble, in- 
cludes the standard statistical ensembles as particular 
cases. 

Let us illustrate above statements for a simple sys- 
tem of non-interacting massless particles, neglecting 
the effects of quantum statistics (Boltzmann approx- 
imation). For a particular realization of Eq. ([1]) we 
choose: 



PaiN) = JdE P^{E) P„ 



(2) 



to calculate the multiplicity distribution Pa{N) in 
the presence of an energy distribution Pa{E). We 
will firstly discuss the MCE multiplicity distribution 
Pmce{N; E) and then solve the integral ^ for a par- 
ticular choice of Pa{E) in the large volume limit. 
The MCE multiplicity distribution is given by 0], 



Pmce{N; E) 



1 
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where x = gVE^/ir^, and g is the particle's 
degeneracy factor. The MCE partition function 
Zmce{y,E) is defined by the normalization condi- 
tion, Y.N=iP^nce{N-E) = 1. In the TL, x » 
1 in Eq.([3]), one finds from the raw moments 
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Additionally it follows from the equivalence of statis- 
tical ensemble that GCE and MCE values for average 
multiplicity are equal to each other, {N)mce = N = 
gVT^/n^, where the temperature T of the GCE is 
found from E = E = igVT^ / . Momentum spec- 
tra in the MCE also converge to GCE Boltzmann 
spectra under this limit. Multiplicity fluctuations, ex- 
pressed by the scaled variances, are however different, 

Ugce = 1 and Wmce = 1/4. 



In the TL the MCE distribution ([3]) converges to a 
Gaussian, 



P„,ce{N;E) 
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The Normal form of distributions is a general feature 
of all statistical ensembles in the TL [11[. 



The GCE energy distribution is equal to: 
1 



PgceiE) 



Zgce{V,T) 



exp 
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where the GCE partition function Zgce{V,T) 
is defined by the normalization condition, 
/ dE Pgce{E) — 1. We consider a system 
with fixed large volume V. As an illustrative example 
we will use the following asymptotic form of energy 
distribution. 



Pa{E) = 



1 
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exp 
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where E — SgVT^/n^ is the GCE energy expectation 
value and uje — {E"^ — E )/E = AT is the scaled 
variance of GCE energy fluctuations. This choice for 
Pa (E) results in a simple correspondence to the GCE 
and MCE in the large volume limit. In Eq. fT]), a is 
a dimensionless tuneable parameter for the width of 
the distribution. In the MCE hmit a ^ 0, Eq. ^ 
becomes a Dirac (5-function, S{E — E). For a = 1, 
Eq. dll) results in the GCE energy fluctuations ([6]) in 
the TL jll|]. The physical interpretation of the GCE 
energy fluctuations corresponds to an 'infinite heat 
bath'. The values, < a < 1, would correspond to 
a 'large, but finite' heat bath. The case of a > 1 we 
would like to denote as 'strong' energy fluctuations. 

Lastly we note that the MCE distribution 
Pmce{N; E) can be conveniently presented in terms 



m, PmceiN;E) 



of the GCE distributions 

Pgce{N,E)/PgceiE). At fixed volume V, the GCE 
distribution Pgce{E) is given by Eq. ([7]) with a = 1, 
while the joint GCE energy and multiplicity distribu- 
(4) tion Pgce {N, E) is given by a bivariate normal distri- 



bution in the large volume limit 11| . 



Pgce{N,E) 
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where AE; = S - AN ^ N - N , E = VKf, N 



nfuJE, 



'N 
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and the correlation coefficient S = (ten / (cte^'n), with 
aEN = ^2'^ ■ The relevant cumulants k, for both Att- 
integrated particle yields and for yields in limited mo- 
mentum windows Ap for the ultra-relativistic Boltz- 
mann gas are given in the Appendix. One finds then 
for the multiplicity distribution: 



1 



PgceiN,E) 

exp 



gee 



, (9) 
(10) 



= 1 and LOmce = 1/4. As it 
for a = 0, and LOa — ujqce 



(27ru;„iV)i/2 

where for our system tOgce 
can be expected, Wq. = bJmc 
for a = 1. It also follows, ujmce < tUa < ^gce for 
< a < 1, and cua > ^gce for a > 1. The a-ensemble 
defined by Eqs. (|2|7p presents an extension of the GCE 
(a — \) and MCE (a = 0) to a more general energy 
distribution. Different values of a correspond, by con- 
struction, to the same expectation values of energy E, 
and multiplicity N . Hence in the TL all ensembles 
defined by Eq.([7| are thermodynamically equivalent. 
Energy and multiplicity fluctuations are however dif- 
ferent. 

As a next step we want to discuss the ef- 
fect of energy fluctuations, Eq.(l7]), on multiplic- 
ity fluctuations in limited segments of momentum 
space for a classical ultra-relativistic gas in the 



TL. The results for the scaled variances, 



{{^Ap)a - {NAp)a) / {N Ap) a , indifferent momentum 
bins Ap are shown in Fig. 1 for several values of 
parameter a. As in Ref.[12|, each momentum bin 
Ap = [pi,P2] contains the same fraction q of the to- 
tal average multiplicity, q = {N^p)a/{N)a. In the 
MCE (a = 0), multiplicity fluctuations are essentially 
suppressed with respect to the GCE {a = 1). Global 
energy conservation in the MCE introduces correla- 
tions in momentum space |l^. The larger the frac- 
tion of the total energy in a given momentum bin, 
the stronger is the MCE suppression effect (a = in 
Fig. 1). For < a < 1 the suppression effects be- 
come weaker, a — 1 corresponds to the GCE results 
for multiplicity fluctuations: uja — I, both in full mo- 
mentum space, and in different momentum bins. For 
'strong' energy fluctuations, a > 1, we find an in- 
crease of multiplicity fluctuations with increasing bin 
momentum, as shown in Fig. 1. 

In order to make a quantitative comparison with 
present and future data on multiplicity fluctuations 
in relativistic nucleus-nucleus collisions one needs to 
extend essentially the formulation considered in this 
letter. Please note that Eqs. (|2l7l l8|) can be readily gen- 
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FIG. 1: Momentum dependence of the scaled variance tj^^ 
for a classical ultra-relativistic gas at T = 160 MeV. Mo- 
mentum bins are constructed in a way that each bin con- 
tains the same fraction g = 0.1 (upper panel), or q — 0.2 
(lower panel) of the average particle yield. The horizontal 
bars indicate the width of bins, while the marker indicates 
the center of gravity of the corresponding bin. Calcula- 
tions are done for different values of a. The total par- 
ticle yield would correspond to g = 1, and according to 
Eq. (|10|l the scaled variance Ua for the total yield fluctu- 
ations equals to 1/4, 7/16, 1, 31/16 at q=0, 0.5, 1, 1.5, 
respectively. 



eralized to more complicated cases of non-zero parti- 
cle masses, quantum statistics, many particle species, 
several conserved charges, and multiplicity fluctua- 
tions in limited segments of momentum space. In 
all these cases the expressions for Pmce{N;A) and 

The MCE 



Pgcei^) are already obtained [ll|, |l 



multiplicity distributions in systems with full hadron- 
resonance spectrum and quantum statistics effects 
were presented in Ref. [Hi] . The role of momentum 
conservation was discussed in Ref. [l^ • The next step 
would be an extension of this formulation to an exter- 
nal distributions PaiA) = P^iV, E, P, B, S,Q), e.g., 
taking into account fluctuations of energy E, momen- 
tum P, and three Abelian charges - baryon number 
B, strangeness S, and electric charge Q. This is to 
be done according to Eq. ([1]). Fluctuations and cor- 
relations of the extensive quantities presented by the 
distribution Pa{V, E, P, B, S,Q) could then be con- 
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nected to measured fluctuations of hadron multiplic- 
ities in limited segments of momentum space. This 
will be a subject of future studies. 

In this letter we suggested to extend the concepts 
of statistical ensembles and introduced the a-ensemble 
defined by an external distribution of extensive quan- 
tities. The key assumption used for calculation of mul- 
tiplicity distributions is the 'equiprobability' of all mi- 
crostates with the same set of extensive quantities. We 
have discussed a simple example of a gas composed of 
classical massless particles to demonstrate the effects 
of energy fluctuations on multiplicity fluctuations in 
full momentum space and in limited segments of mo- 
mentum space. Measurement of event-by-event multi- 
plicity fluctuations in different momentum bins corre- 
spond to current experimental studies of hadron pro- 
duction in relativistic nucleus- nucleus collisions. We 
believe also that the concept of statistical ensembles 
with fluctuating extensive quantities may be appropri- 
ate in other situations too. In fact, in all cases when 
fluctuations of extensive quantities are a subject of 
interest and can be measured experimentally. 
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Appendix 

For the ideal Boltzmann gas of massless particles 
the momentum distribution is f{p) = cxp(— p/T), and 
the cumulants are: 



E,N 
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«f = dpp' fip) 



gT^ 



= 3 



Additionally, in Boltzmann approximation, ' 

and K^'^ = k]^. The cumulants in the momentum 
segments Ap — [p2,Pi] are: 
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Additionally, in the Boltzmann approximation 
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